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Abstract 

In this paper we examine the Casimir effect for charged fields in presence of 
external magnetic field. We consider scalar field (connected with spinless parti- 
cles) and the Dirac field (connected with 1/2-spin particles). In both cases we 
describe quantum field using the canonical formalism. We obtain vacuum energy 
by direct solving field equations and using the mode summation method. In order 
to compute the renormalized vacuum energy we use the Abel-Plana formula. 



1 Introduction 

The imposition of boundary condition on a quantum field leads to the modification of 
the vacuum energy level and can be observed as an associated vacuum pressure. This 
effect (called Casimir effect) has been predicted by Casimir in his original work in 
1948[1], and has been experimentally observed for electromagnetic field several years 
later. Until now, in many theoretical works Casimir energy has been computed for 
various types of boundary geometry and for fields other than electromagnetic one. 

When we consider of charged fields another important question arise. How an 
external field, coupled to the charge, affect the vacuum energy of the field? Answer 
to this question is important for understanding of some aspects in particle physics. 
Within a hadron, for example, the vacuum energy of quark fields is affected by the 
electromagnetic field of the quarks and by the color field of gluons and quarks. 

In order to investigate how charged fermionic and bosonic constrained vacuum 
fluctuations are affected by fields coupled to this charge, we consider vacuum energy of 
electrically charged fields under the influence of an external constant uniform magnetic 
field, and constrained by simple boundary conditions. In Section 2. we shall consider 
a complex scalar field confined between two infinite plates with Dirichlet boundary 
conditions with magnetic field in a direction perpendicular to the plates. In Section 3. 
we shall consider a Dirac quantum field under antiperiodic boundary conditions. This 
choice of geometry and external fields avoids technical difficulties and focuses our 
attention on the fundamental issue. 

The fermionic Casimir effect was first calculated by Johnson [2] for applications 
in the MIT bag model. For a massless Dirac field, Johnson obtained an energy density 
7/4 times the energy density of the electromagnetic field. 
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In both cases (bosonic and fermionic) we describe quantum field using the canon- 
ical formalism. We obtain the vacuum energy by direct solving the field equations 
and using the mode summation method. In order to compute the renormalized vacuum 
energy we use the Abel-Plana formula. 

In order to obtain influence of magnetic field on the Casimir energy another meth- 
ods can be used. Recently, the papers using Schwinger's proper time method have been 
published[3, 4]. So, our purpose is to compare results obtained using our method with 
the Schwinger's one. 

2 Bosonic Casimir effect 

In this section we consider scalar field 0(r, t), describing charged, spinless particles 
with mass m. It is under the influence of uniform magnetic field B = (0,0, B). We 
choose direction of z axis in such a way that B is positive. 

Using gauge invariance, we choose electromagnetic potential corresponding to field 
Bin the form: A M = (0,-yB,0,0). 

The lagrangian density of field takes the form: 

i? = (D M 0)*(Z)^)-mV4> = 

= d t <j)*d t <j)-(d x + ieBy) <j)*(d x - ieBy) <j> + 

-d y <j)*dy<j)-d z <j)*d z <j>-m z <j)*<j), (1) 

and leads to the following equation of motion 

d 2 (p - d 2 <p - dfo - d 2 <j> + 2ieByd x <j> + (e 2 B 2 y 2 + m 2 ) <p=0. (2) 

where derivative D^ = + ieA^ in Eq. (1) is covariant derivative. 

Solution of Eq. (2) is well known and is given in [6] in case of nonrelativistic 
quantum particle. We briefly remind method of solving this kind of equation. The 
variables x, z and t do not occur in Eq. (2) explicity, therefore the solution of Eq. (2) 
takes the form: 

<t>{x,y,z,t)=F{y)e ii t x+pz - m \ (3) 
After inserting of Eq. (3) into Eq. (2) the equation for function F(y) is given by: 

d 2 F + (CO 2 -k 2 -p 2 - m 2 + 2eBky - e 2 B 2 y 2 )F = 0, (4) 

which has a solution in the form: 

F(y)^e- eB / 2 ^ 2 U n (V^B(y-y Q )), (5) 

where H„ are Hermite polynomials and yo = k(eB)~ l . The parameters n, p i CO fulfil 
following relation 

co = ^2eB{n + j)+p 2 + m 2 . (6) 
Let us introduce the Dirichlet boundary conditions in the form: 

<j)(x,y,z = 0,t) = <l)(x,y,z = a,t)=0, (7) 
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which implies that field is equal to zero on the two parallel planes (plates). The plates 
are perpendicular to the z-axis and the distance between them is equal to a. 
In the area between plates, solution of Eq. (2) is given by: 



0(r,f) = r £ (a n ,kUnik(7)e-"°^ + b+ k u nlk (r)e i<0 -") . 

J -°° /=ln=0 



(8) 



Functions u„i k are special solutions of the Eq. (2) in the form: 

unikif) = C n e ikx sm(p lZ )e- eB ^y-y^}i n (V^B(y-y )), (9) 

where C„ = (eBa 2 ) l / 4 (2"n\7t 3 / 2 )~ 1 / 2 is normalisation constant, and pi — nl/a. Solu- 
tions of Eq. (9) fulfil relations of orthogonality in the form: 

J^d 3 ru ni i lkl (r)u* n2l2k2 (r) = a 2 d(ki - k 2 )8i l i 2 8 nin2 (10) 

and relations of completeness in the form: 

/oo co oo 

<**L L Unik(nK lk {r 2 ) =a 2 5 3 (n -r 2 ), (11) 
00 /=ln=0 

where T is area between the plates. 

Standard commutation relations for the field 

Wi,*)M + (h,t)] = i8 3 (7i-?2) (12) 
[0(n,f),0(r 2 ,f)] = [^0+(n,f),^0+(r 2 ,f)]=O (13) 

lead to the following commutators for creation and annihilation operators 

[a ni h{h),a+ 2h {k 2 )] = [b nih {h),b+ 2h {k 2 )} = \a- 2 G>-^8 nin2 8 hh 8{ki-k 2 ). (14) 
The hamiltonian density: 

J? = a t 0*a t + (d x + ieBy)<$>*{d x - ;eBy)0 + <9 y 0*^0 + d z 0*d z +m 2 0*0 (15) 
leads to the vacuum energy in the form: 

E vac =L 2 dz(0\J?\0)=L 2 eB(2ny l Y J Y<\Ji l2 + eB ( 2n + 1 )+ m2 > ( 16 > 

^ 1=1 n=0 " 

where L is length of the plates. 

The energy (16) is infinite. Now, we define the renormalised vacuum energy E ren 
as a difference between energy (16) and energy without boundary conditions: 

E ren = L 2 eB(2n)- 1 £ f£ ^l 2 +eB(2n + 1) +m 2 + 

n=0 V=l 

-jf c//y / ^/ 2 +eB(2n + l)+m 2 ^. (17) 
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In order to compute the renormalised energy we use the Abel-Plana formula in the 
form 

if (I) - [f{l)dl = 4/(0) + i [ 'dt f W- f _ { - U) . (18) 

This formula is often applied in renormalization problems. Its applications and gener- 
alizations can be found in many works [7, 8]. 

The renormalised energy, after using formula (18), takes the final form: 

E ren = -L 2 eBa- x £ f dt ^ t 2 - X 2 (e 2m - l) (19) 

11=0 •> ^" 

where X„ — a%~ 1 ^ eB(2n + 1) + m 2 . 

Numerical calculations of Eq. (19) are shown in Figs. 1-5. Main conclusion is that, 
external magnetic field decrease the Casimir energy and pressure of the scalar field 
(Fig. 1), and suppress it completely in the limit B — > 0. The decrease of the Casimir 
energy with a distance is faster for stronger magnetic field (Fig. 2). 

Let us to compare the behaviour of the Casimir pressure for two different distances 
a (Figs. 3,4). In order to suppress the pressure noticeably, magnetic field should be of 
order of 10r for distance a = 10«m, and of order of 0. 1 T for distance a — lOOnm. 

Asymptotic behaviour of energy (19) for strong magnetic field (B — > °o) we obtain 
using the relation exp(2^f) — 1 — ► exp(27Cf ), which gives: 

E B ^ = -L 2 eBa- 1 £ / dtJt 2 - X 2 e- 2jlt = 

= -\L 2 eB%- 2 £ ^m 2 + eB(2n+l) Kj (2a ^m 2 + 2eB(2n + 1 )) , (20) 

17=0 

where (x) is Bessel function. 
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Fig. 1. Casimir energy for scalar field as a function of magnetic field eB, corresponding to 
a = 1 and L = 1. Bold curve corresponds to mass m = whereas thin curve corresponds to mass 
m=\. The dashed line corresponds to value — jt 2 /720. (It is standard value of energy for B = 0) 
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Fig. 2. Casimir energy for scalar field as a function of distance a, for massles field m — 
and L = 1 . Bold curve corresponds to magnetic field eB = 8 whereas thin corresponds to eB=\. 
The dashed curve corresponds to function — ?r 2 a~ 3 /720 (standard function of distance a forB = 0). 
Values of energy are shown in logarithmic scale. 




Fig. 3. Casimir pressure for scalar field as a function of magnetic field B, corresponding to 
distance a= lOnm. The diagram is made in SI units for massless field. 




Fig. 4. Casimir pressure for scalar field as a function of magnetic field B, corresponding to 
distance a= lOOnw. The diagram is made in SI units for massless field. 
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Fig. 5. Casimir pressure for massless, scalar field as a function of distance a. Diagram is 
made in SI units, using logarithmic scales. Bold curve corresponds to B = 107\ thin to B = \T, and 
dashed to B = 0T. 

Asymptotic behaviour of expression (19) for weak magnetic field (B — ► 0) we ob- 
tain, when we change sum into integration, which gives: 



Eb^o = -L 2 eBcT l 



dn 



X{n) 



dtJt 2 -X 2 (n)(e 2nt -l) = 



-L 2 an- 2 dyy/y 2 - 1 / dxx 3 (e 2axy - 1 ) 



M 



(21) 



where M=y/m 2 + eB. ForM = 0Eq. (21)takesthe standard value £(0) = -L 2 % 2 j (720a 3 ) 
for charged, scalar field [7]. 

We compared Eq. (19) with Eq. (15) from [3] (based on Schwinger's method) using 
Mathematica 4.1. We obtained the same numerical results. However, we are not able 
analytically prove that these two equations are equivalent. 
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3 Fermionic Casimir effect 

The lagrangian density of the Dirac field in external electromagnetic field takes the 
form: 

JSf = ^ + f(ifd^-efA^-m)^, (22) 
and leads to the following equation of motion 

(i'fdp-e'fA ll -m)W = 0, 

where a = 1,2,3,4 are numbers of bispinor components. 
We take y matrices in Dirac representation: 



(23) 



I 

-/ 



o & 

-a 1 



(24) 



Analogously as in previous section, we solve Eq. (23) in presence of an external 
magnetic field B = (0,0, B), but now we impose antiperiodic boundary conditions in 
the form: 

V(x,y,z = 0,?) = -V(x,y,z = a,t). (25) 



Solution of Eq. (23) can be predicted as follows: 

t) = exp(i(fiM + kx + pz))ri a (y), 
where bispinors T] a (y) take the form: 



(26) 



ri a (y) 



(ciu n -i(y)\ 
ciu n {y) 

\ c 4 u„(y) J 



(27) 



where c a are number parameters independent of space co-ordinates, and functions 
u n {y) are given by: 



u n {y) =e- eB Wy-y^H n (V^B(y~y )), 

where yo = k/(eB). 

After inserting Eq. (28) into Eq. (23) we obtain matrix equation: 



(28) 



/-(co + m) 





~P 


2nVeB\ 




f Cl \ 





-(co+m) 


VeB 


P 




C2 


P 


-2nVeB 


co — m 









\ -\feB 


-P 





co — m) 




w 



(29) 



Nontrivial solutions of Eq. (29) (i.e. different from null vector) exist if rank of matrix 
is different than 4. This situation occurs when: 



CO 2 = p 2 + 2eBn + m 2 . 



(30) 
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After inserting solutions of Eq. (29) into Eq. (26) we obtain general solution of 
Eq. (23) in following form 



poo , _^ 

W = EE / ^E(^ ,(£r+fa+m) K./,«^./,(^y)+ 

n =Qi=-o° J ~°° .v v 

+e i{Et+ ^ +Plz) )b+ l /k)n n ,iAk,y)), (3D 



where /?/ = n/a(2l + 1), £ = |o|, and variable s takes the values -1, 1. 
For n > bispinors 77 from Eq. (31) are given by: 



ln,s=-i{k,p,y) = 



r\n, s =i{k,p,y) 



nn,s=-i{k,p,y) = 



(--^ u n-i(y)\ 

u„-i(y) 




r 



/2"n! 



f 2 -w^»n-i(y)\ 



E-\-m 



V un(y) J 

( 

u„(y,-k) 

■ 2j ^un-i{y,-k) 

V 7^E u n(y- k ) I 



^/2"-\n-\)\ 



( u n _ x {y,-k) \ 


-^E U n-\{y-k) 



\ 



(32) 



(33) 



(34) 



(35) 



eB u n (y,-k) J 



m+E 



where C„ = ((eB / 'n) l l 2 {E + m) /2E) 1/2 . 

For n = bispinors <^o,.s=-i an d ^)o,j=i are a l so given by these equations, if we take 
M -i (y) — 0- Whereas spinors 77o,.s=-b <^o..s=i are equal to zero. 

The spinors 77, E, are orthonormal in the sense of the norm defined by: 

||^(*,p,y)|| 2 = dytAk,P,y)&, s (k,p,y). (36) 

./ — CO 

Anticommutation relations for the field operators X P, II = fV + in the form: 

\V a (x,t),Il b (y,t)} + = i8(x-y)8 ab (37) 
lead to the following anticommutators for creation and annihilation operators 

[a ni h Sl {ki),a^ hS2 {k 2 )) + = [b ni ^ Sl {k\ ) , b„ 2l2S2 {ki)]+ = 

= (2^/nay l 8 S[S2 8i [ i 2 8 nin2 8{ki -k 2 ). (38) 
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In Dirac theory, the vacuum state is defined as the state with filled up negative 
energy levels. This situation corresponds to: 

a nls \0) =0 b nls \0) =0. (39) 

The hamiltonian density of the field: 

St? = -*i> + f(iYdi - efA^ - mfP = f¥+%V (40) 

(where in second step we use the fact, that VP fulfils the Dirac equation) leads to the 
vacuum energy density as follows 

S vac (a) - (0| J^O) = 

= ijk L El/ ^£(4VMWM+n«Vi(Mn«/,-i(M) = 

n=0/=-=o s - 7 - 00 



= -eB{y/na)- 1 £ £ a„y ^-(2/ + l) 2 + 2eBn + m 2 (41) 

„=0/=-oo 

where OCo = 1/2 and CC„ = 1 for n > 0. 

The energy density (41) is infinite. Renormalised vacuum energy is difference 
between energy density (41) and vacuum energy without boundary conditions, and is 
given by: 

S ren (a) = -2eB{Vna)- 1 £ Cfe(£ J \% 2 a - 2 {l + l/2) 2 + 2eBn + m 2 + 

n=0 \=0 

-J dl \/\% 2 a- 2 l 2 + 2eBn + m 2 j . (42) 
In order to calculate Eq. (42) we use Abel-Plana formula in the usefull form: 

t (l - I 12)- I (III (I) = -i I — "' 

t=0 Jo Jo e 1 

and finally we obtain renormalised energy density in the form: 



^F{t + \/2)- T dtF{t) = -i T ^—(F(it)-F{-it)), 
r~k Jo Jo e Lm + 1 



(43) 



,(«) = SeBa- 2 £ a n f V^^' 

n=0 JA " c 1 



(44) 



where X n — a(2n)~ x ^f2eBn + m 2 . 

Numerical calculations of Eq. (44) are shown on Figs. 6,7. In contrast to the 
bosonic case, for the Dirac field the Casimir energy is enhanced by the external mag- 
netic field (Fig. 6). So, the Dirac vacuum behaves like paramagnetic medium. Decrease 
of the Casimir energy with a distance is slower in strong magnetic field (Fig. 7). 
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Fig. 6. Casimir energy for the Dirac field as a function of magnetic field eB, corresponding to 
distance a = 1 and L = 1. Bold curve corresponds to mass m = whereas thin corresponds to mass 
m = 1. Dashed line corresponds to value equal to -7jt 2 /180 (standard value of energy for B = 0). 
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Fig. 7. Casimir energy for the Dirac field as a function of distance a, for massles field m = 
and L= 1. Bold curve corresponds to magnetic field eB = 8 whereas thin corresponds to eB = 1. 
Dashed curve corresponds to function — 7?r 2 a~ 3 /180 (standard function of energy for B = 0). Values 
of energy are shown in logarithmic scale. 

Vacuum pressure is shown on Figs. 8-9. For massless field (and strong B) it is 
increasing linear function of magnetic field B (Fig. 8). For massive field, pressure 
depends on the field B very weakly - two curves on Fig. 9 (for B — IT and B = \QQT) 
coincide. For electron field values of Casimir pressure are measurable for distances 
smaller then Q.lnm (Fig. 9). In practice it excludes possibilities of experimental tests. 
It can play a role only in subatomic distances, for example in some particle models. 

Analogously as for bosons, we compared Eq. (44) with Eq. (9) from [4] (based on 
Schwinger's method) using Mathematica 4.1. We obtained the same numerical results. 
However, we are not able analytically prove that these two equations are equivalent. 
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Fig. 8. Casimir pressure for massles, the Dirac field as a function of magnetic field B, 
corresponding to distance a = lOOnm. The diagram is made in SI units. 




Fig. 9. Casimir pressure for the Dirac field of electrons (m = m e ) as a function of distance a. 
The diagram is made in SI units for two values of field B= IT and B = 1007\ Both curves coincide. 

Asymptotic behaviour of energy density (44) for strong magnetic field (B — > °°) 
and for am — > °° we obtain using the relation exp(27£f) — 1 — > exp(27£f), what gives: 

S B ^ = -8eBa- 2 £ a„ f dt J ' t 2 - X 2 e' 27 " = 

= -leBn^cT 1 «« V m 2 + 2eBnKi (ay/m 2 + 2eBn) . (45) 

n=0 

Asymptotic behaviour of Eq. (44) for weak magnetic field (B — > 0) we obtain 
changing sum into integration, which gives: 

g B ^ Q = -%eBa- 2 J dnj dt ^t 2 - X 2 (n) (e 2m + l) = 

= -2%- 2 l~ dyy/y 2 -\ (°° dxx 3 {e 2ax7 + l) -1 , (46) 

J l Jm 

For m = expression (46) is equal to <^(0) = — 77T 2 /(180a 4 ), which is standard value 
of Casimir energy of the Dirac field without external fields [2]. The expression (46) 
donot depends on B, which suggests that Ss-tQ = <^(b=0) + 0{B 2 )- 



11 



Acknowledgements 



I wish to thank B. Broda for numerous discussions about results. 
This article was supported by the Polish Ministry of Scientific Research and Informa- 
tion Technology under the grant No. PBZ/MIN/008/P03/2003 and by the University of 
Lodz. 



References 

[1] H. B. G. Casimir, Proc. Ron. Nederl. Akad. Wet. 51, 793 (1948). 
[2] K. Johnson, Acta Phys. Polonica B6 865 (1975). 

[3] M. V. Cougo-Pinto, C. Farina, M. R. Negrao and A. C. Tort, "Bosonic Casimir 
effect in external magnetic field", J.Phys. A32, 4457-4462, (1999). 

[4] M. V. Cougo-Pinto, C. Farina and A. C. Tort, "Fermionic Casimir effect in an 
external magnetic field", Braz. J. Phys. 31 84 (2001). 

[5] E. Elizalde, F. C. Santos and A. C. Tort, J. Phys. A 35 7403 (2002). 

[6] L. Landau and E. Lifszic, Quantum Mechanics: Non-Relativistic Theory, (Perga- 
mon, Oxford, 1977). 

[7] M. Bordag, U. Mohideen and V. M. Mostepanenko, "New developments in the 
Casimir effect", Phys. Rept. 353 1-205 (2001). 

[8] A. A. Saharian, "The generalized Abel-Plana formula. Applications to Bessel 
functions and Casimir effect", e-print hep-th/0002239 (2000). 

[9] A. A. Grib, S. G. Mamayev and V. M. Mostepanenko, Vacuum Quantum Effects 
in Strong Fields, (Friedmann Laboratory Publishing, St. Petersburg, 1994). 

[10] V. M. Mostepanenko and N. N. Trunov, The Casimir Effect and its Applications, 
(Clarendon Press, Oxford, 1997). 

[11] L. H. Ford, Phys. Rev. D 21 933. 

[12] M.A. Johnson and B.A. Lippmann, Phys. Rev. 76 , 282 (1949). 



12 



